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Transverse-momentum resummation for gaugino-pair production at hadron colliders 
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We present a first precision analysis of the transverse-momentum spectrum of gaugino pairs pro- 
duced at the Tevatron and the LHC with center-of-mass energies of 1.96 and 10 or 14 TeV, respec- 
tively. Our calculation is based on a universal resummation formalism at next-to-leading logarithmic 
accuracy, which is consistently matched to the perturbative prediction at 0(a s ). Numerical results 
are given for the "gold-plated" associated production of neutralinos and charginos decaying into three 
charged leptons with missing transverse energy as well as for the pair production of neutralinos and 
charginos at two typical benchmark points in the constrained MSSM. We show that the matched 
resummation results differ considerably from the Monte Carlo predictions employed traditionally in 
experimental analyses and discuss the impact on the determination of SUSY mass parameters from 
derived transverse-mass spectra. We also investigate in detail the theoretical uncertainties coming 
from scale and parton-density function variations and non-perturbative effects. 



INTRODUCTION 

Weak-scale supersymmetry (SUSY) is a very well mo- 
tivated extension of the Standard Model (SM) of parti- 
cle physics. It can break the electroweak symmetry ra- 
diatively, allows for its grand unification with the local 
gauge symmetry of strong interactions, offers a natural 
explanation of the large hierarchy between electroweak 
and gravitational interactions, and appears naturally in 
string theories [1]. Among the new particles predicted by 
the Minimal Supersymmetric SM (MSSM), the fermionic 
partners of the neutral and charged gauge and Higgs 
bosons, called neutralinos and charginos, may be rel- 
atively light, and the lightest neutralino, stabilized by 
an at least approximate i?-symmetry, represents one of 
the most promising dark matter candidates, whose gaug- 
ino/higgsino decomposition has important consequences 
for cosmology. The search for SUSY particles and the 
identification of their properties have thus become defin- 
ing tasks of the current hadron collider program. Par- 
ticular attention has since long been paid to the produc- 
tion of gauginos [2], which are produced either directly 
or through squark/gluino decays, decay themselves le 
tonically and are easily identifiable at the Tevatron 
and at the LHC dj. The Tevatron collaborations CDF 
and DO have already published several SUSY searches 
in the "gold-plated" trilepton channel For the LHC, 
the CMS collaboration estimate the reach of the direct 
channel to be relatively modest (universal gaugino mass 
parameter m 1 / 2 < 180 GeV) and relevant only for scenar- 
ios with heavy colored particles (universal scalar masses 
mo > 1000 GeV). However, the ATLAS collaboration 
have shown that a discovery can be made in this channel 
for a wider range of gaugino masses {m 1 / 2 < 360 GeV) 
and independently of mo already with a low luminosity 
of a few fb _1 [4J. Gaugino pair production is therefore 



a very important SUSY discovery channel at both cur- 
rently running hadron colliders. 

For an efficient suppression of the SM background from 
vector-boson and top-quark production and a precise de- 
termination of the underlying SUSY-breaking model and 
masses, an accurate theoretical calculation of the signal 
(and background) cross section is imperative. As the 
lightest SUSY particle (LSP) escapes undetected, the key 
distribution for SUSY discovery and measurements is the 
missing transverse-energy ($t) spectrum, which is typi- 
cally restricted by a cut of 20 GeV at the Tevatron and 
30 GeV at the LHC. While the SUSY particle pair is pro- 
duced with zero transverse momentum {pr) in the Born 
approximation, the possible radiation of gluons from the 
quark-antiquark initial state or the splitting of gluons 
into quark-antiquark pairs at 0(a s ) in the strong cou- 
pling constant induces transverse momenta extending to 
quite substantial values and must therefore be taken into 
account In addition, the perturbative calculation di- 
verges at small px, indicating the need for a resummation 
of soft-gluon radiation to all orders @ . Only after a con- 
sistent matching of the perturbative and resummed cal- 
culations an accurate description of the (missing) trans- 
verse energy spectrum and precise measurements of the 
SUSY particle masses can be achieved. 

In this Letter, we report on the first precision analysis 
of the transverse-momentum spectrum of gaugino pairs 
produced at the Tevatron and the LHC with center-of- 
mass energies of 1.96 and 10 or 14 TeV, respectively. We 
briefly describe in the following section our implemen- 
tation of the resummation formalism, which has been 
improved with respect to the original proposal in nu- 
merous respects, and present then numerical results for 
the production of various gaugino pairs at two typical 
MSSM benchmark points. We also discuss the impact 
of the computed precise transverse-momentum spectrum 
on the determination of SUSY mass parameters and in- 
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vestigate in detail the remaining theoretical uncertainties 
coming from scale and parton-density function variations 
and non-perturbative effects. 



TRANSVERSE-MOMENTUM RESUMMATION 

In the Born approximation, the production of neutrali- 
nos and charginos at hadron colliders 



PP, PP^ <?<?' + x -> Xi X] + X 



(1) 



is induced by the quarks q and antiquarks q' in the initial 
(anti-)protons and is mediated by s-channel electroweak 
gauge-boson and t- and it-channel squark exchanges. Its 
partonic cross section a^ ' can be expressed in terms of 
the gaugino and squark masses to^o,± and rriq, the masses 
of the electroweak gauge bosons, the Mandelstam vari- 
ables s, t and u, and generalized charges @]. 

At leading order (LO) in the strong coupling constant, 
0{a s ), virtual loop and real parton emission corrections 
must be taken into account [5[. The latter induce trans- 
verse momenta of the gaugino pair, that extend typically 
to values of the order of the gaugino mass. In the small- 
Pt region, where the bulk of the events is produced, the 
convergence of the perturbative expansion is spoiled due 
to the presence of large logarithms a"/p^ \n m (M 2 /p T ) 
with m < 2n — 1. These must be resummed to all orders 
in impact parameter (b) space, 



da 



RES 



dp. 



<PT,M, s) 



M 2 



db-J (b PT )W(b,M, s ), 



(2) 

in order to correctly implement transverse-momentum 
conservation. Here, b describes the minimal distance of 
the two incident particles in the limit of no interaction 
and is the conjugate variable of the transverse momen- 
tum pt; Jo(x) is the O'^-order Bessel function; M is the 
invariant mass of the gaugino pair; and Eq. (2) is evalu- 
ated numerically as in Ref. [7] . After a Mellin transform 
with respect to z = M 2 /s, the hadronic cross section 
simplifies from a convolution to a product of the parton 
densities evaluated at the factorization scale [ip with the 
partonic cross section, the TV-moment of the latter 

W N (b, M) = U N (M) exp [Qn(L)] (3) 

factorizing further into a 6-independent function 



H N {M) = a w (M) 1 + — H%\m/Q) + 

that is therefore finite as pr — > or b — > 
exponential form factor 

Qn{L) = LgW(a s L) + g (2 \a s L) + 



(4) 



oo, and an 



(5) 



which resums the divergent leading and next-to-leading 
contributions in the logarithm L — \n(Q 2 b 2 /6q) through 



the process-independent functions g^ 1 ' 2 ^ ■ While &o = 
2 e -7-B j s f kinematical origin, the scale Q ~ M accounts 
for the arbitrary separation of the two factors. Unphys- 
ical logarithmic divergences at b — > are regularized by 
replacing L with L = \n{Q 2 b 2 jb\ + 1). The evolution 
from the factorization scale hf and the low scale bo/b to 
the high scale Q is performed in Hn(M) and exp[(?jv], re- 
spectively, leaving the latter [ip -independent and allow- 
ing for a convolution of the partonic cross section with 
the parton densities in the (anti-)protons at p,p 

Although SUSY particles appearing in the virtual cor- 
rections must be taken into account for a proper renor- 
malization of ultraviolet divergences, their heavy masses 
are known to leave only little imprint numerically. In par- 
ticular, the ratios of NLO/LO total cross sections vary 
only from 1.26-1.28 at the Tevatron and 1.24-1.28 at the 
LHC for squark masses of 350-1000 GeV (for the similar 
case of sleptons, see also Fig. 4(a) in Beenakker, Klasen 
ct al. [5]). The restoration of the equality of the weak 
scalar (<?) and vector gauge couplings (g) furthermore re- 
quires the introduction of a finite SUSY-restoring counter 
term, g = g[l — a s /(67r)]. Finally, quark-gluon initial 
states can induce the production of gauginos in associa- 
tion with real squarks decaying into quarks and gauginos. 
The singularity associated with the pole of the on-shell 
squark propagator is regularized by a small finite width 
Tq ~ 10~ 2 m,j (the exact value has little influence numer- 
ically), and the on-shell contribution is subtracted from 
the total cross section to avoid double-counting 0, Q . 



NUMERICAL RESULTS 

To obtain a valid hadronic cross section at all values of 
Pt, the 0(a s ) (LO) and resummed (RES) partonic cross 
sections are matched by subtracting from their sum the 
perturbatively expanded (EXP) resummed cross section, 



da 

dp 2 



da 



LO 



dpi 



da 



RES 



da 



EXP 



dpi 



dpi 



(6) 



and by performing numerically the necessary inverse 
Mellin and Fourier transforms, kinematic integrations, 
and parton density convolutions. The parton densities 
are evaluated in the most recent parameterization by 
the CTEQ collaboration CTEQ6.6M [§] with ^ F (and the 
renormalization scale /ir) set to the average mass to of 
the final state particles. The SUSY particle masses are 
obtained from universal parameters defined at the grand 
unification scale through the renormalization group run- 
ning implemented in the computer code SUSPECT2 . 3 [l(| ■ 
To be specific, we choose two minimal supergravity 
benchmark points. The first is the low-mass point LMO 
(SU4) with universal fermion mass mi/2 = 160 GeV, 
scalar mass too = 200 GeV, trilinear coupling Aq = —400 
GeV, bilinear Higgs mass parameter /i > 0, and ratio of 
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Higgs vacuum expectation values tan/3 = 10 [4]. It has 
been been denned by the CMS (ATLAS) collaboration 
with the objective of high cross sections and thus early 
discovery at the LHC, as the resulting gaugino, squark 
and slepton masses m^o = m-± = 113 GeV, m^o =61 
GeV, m~ q ~ 420 GeV, and mj ~ 220 GeV lie just beyond 
the current Tevatron limits. In this scenario, the lightest 
chargino and second-lightest neutralino decay with 35% 
and 15% probability through virtual sleptons to the LSP 
and one and two charged leptons, respectively 



pp -> x 



f X ° at = 1-96 TcV - LMO 



11] 



Our second benchmark point is the Snowmass point 
(and slope, SPS) la' with to 1/2 = 250 GeV, m = 70 
GeV, A = -300 GeV, fi > 0, and tan ,3 = 10, which is 
similar to the point LM1 (SU1) defined by the CMS (AT- 
LAS) collaboration and, with its SUSY particle masses of 



m 



m-± 



98 GeV, 



550 GeV, 



and m ; ~ 



183 GeV, to. 
120 . . . 190 GeV, is known to be compatible 
with all high-energy mass bounds and low-energy preci- 
sion data. Here, the lightest chargino and second-lightest 
neutralino decay almost exclusively to three charged lep- 
tons and missing transverse energy, albeit through real 
sleptons, which may be experimentally reconstructed 
through cndpoints in kinematic distributions 12[. 

In Fig. [U we show the corresponding transverse- 
momentum spectra of chargino-neutralino pairs pro- 
duced at run II of the Tevatron (top) and the initial 
run of the LHC (bottom) with center-of-mass energies 
of = 1-96 and 10 TeV, respectively. As expected, the 
LO predictions (dashed curves) diverge at lowpr, but be- 
come finite after matching them to the resummed predic- 
tions at next-to-leading logarithmic (NLL) accuracy (full 
curves). In this region, the perturbative expansions of the 
resummed predictions (dots) coincide with those at LO, 
while at large px they coincide with the resummed ones. 
Through resummation, the perturbative predictions are 
considerably enhanced even at values of px, which are 
of the order of the experimental I$x cuts. It is therefore 
important to clearly distinguish the effects induced by 
QCD radiation and by the unobserved LSPs and neutri- 
nos. By construction, the matched LO+NLL prediction 
allows to reproduce the correct 0(a s ) correction (K) fac- 
tor of the total perturbative cross section after integra- 
tion over pt, e.g. of 1.26 at the Tevatron. For comparison, 
we also show the matched LO+NLL p^-spectrum (dot- 
dashed curve) for the 14 TeV design energy of the LHC, 
which extends to considerably larger values of pt than 
at 10 TeV. The theoretical predictions are influenced by 
three main sources of uncertainty: scale variations, evalu- 
ated in the canonical range of /Uf,.r/to = 0.5 ... 2 (shaded 
bands), variations of the parton densities, evaluated 

through Ap DF = \Jj2i=i(^ a t ~ dcr~) 2 /(2d<7) along the 
22 eigenvector directions defined by the CTEQ collabora- 
tion (lower curves), and three choices of non-perturbative 
(NP) form factors, evaluated through Anp = (dcr NP — 
da) /da (inserts) [13j. For px > 5 GeV, all theoretical 
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FIG. 1: Transverse-momentum spectra of chargino-neutralino 
pairs at the Tevatron (top) and the LHC (bottom). The LO 
calculation (green/dashed) is matched to the resummed cal- 
culation (red/full) by subtracting its fixed-order expansion 
(dotted). The scale uncertainty is shown as a shaded band, 
the PDF (below) and non-perturbative (insert) uncertainties 
as separate graphs, and the matched result for the LHC design 
energy of yS = 14 TeV as a dot-dashed line (bottom). 



uncertainties are smaller than 5% for the LO+NLL pre- 
dictions. In particular, the 5% PDF uncertainty is similar 
to the one obtained for weak boson production 

In Fig. [3J the matched py-spectra for chargino (dotted) 
and neutralino (dot-dashed) pairs are compared to those 
of the tri-lepton channel (full/dashed) discussed above. 
While positive and negative charginos are produced with 
equal rates in pp collisions at the Tevatron, their rates 
differ slightly in pp collisions at the LHC. The cross sec- 
tions for neutralino pair production are about one order 
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FIG. 2: Transverse-momentum spectra at LO+NLL for the 
associated production of charginos and neutralinos (full and 
dashed) as well as chargino (dotted) and neutralino (dot- 
dashed) pairs (dotted) in three different collider modes. 
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FIG. 3: Transverse-momentum spectra at LO+NLL (full), LL 
(dashed), and generated by the PYTHIA parton shower with 
default (bars) and tuned (crosses) parameters at the Tevatron. 



of magnitude smaller, as the second-lightest neutralino 
couples to the s-channel Z°-boson only through its rela- 
tively small Higgsino component. 

In experimental analyses, QCD radiation in hadronic 
collisions is usually simulated with tree-level matrix ele- 
ments and parton showers based on an exponential Su- 
dakov form factor, which resums the leading logarithms 
(LL) and some next-to-leading logarithms. In Fig. 02 
we compare therefore our matched LO+NLL prediction 
(full curve) with our resummed prediction at LL order 



(dashed curve) and the default (bars) and tuned (crosses) 
predictions of the PYTHIA6 . 4 Monte Carlo (MC) gener- 
ator While the default MC prediction is clearly 
improved beyond the LL approximation and approaches 
the LO+NLL result, it peaks at too small values of pr- 
Tuning the intrinsic px of the partons in the proton to 2.1 
GeV for Z°-bosons (CDF tune AW) and 4 GeV for gaug- 
ino pairs (our tune AW) improves the description of the 
peak, but still underestimates the intermediate pr-region 
and the mean value of p T (14 GeV for PYTHIA6 . 4, 15 GeV 
for our tune AW, and 18 GeV for our LO+NLL predic- 
tion). This has, of course, a direct impact on the de- 
termination of the gaugino (and slepton) masses through 
variables derived from the transverse momenta of the ob- 
served leptons pt.i and $t, such as the effective mass 
Msff = J2iPT,i+ $t [3 or the stransverse mass [l6j . 
A detailed discussion of signal and background distribu- 
tions in these variables, including also the experimental 
resolution, is beyond the scope of this work. Let us men- 
tion, however, that the contribution of unmeasured (low- 
Pt or forward) or mismeasured hadronic energy to the 
"fake" $t is under close scrutiny both at the Tevatron 
and at the LHC. The ATLAS trilepton analysis, e.g., does 
not identify jets with px < 10 GeV, and an optional cut 
on jets with px > 20 GeV reduces the significance con- 
siderably. As the two LSPs are often back-to-back, the 
f,T in the trilepton analysis is required to be relatively 
small (> 30 GeV). It can then be affected by an error of 
up to 10% (Aad et al. ||). 

CONCLUSION 

In summary, we have calculated the transverse- 
momentum spectrum of gaugino pairs produced at 
hadron colliders at next-to-leading logarithmic accuracy. 
We have demonstrated that this renders the perturbative 
prediction finite, modifies considerably the traditionally 
used Monte Carlo predictions, and reduces the remain- 
ing theoretical uncertainties to the level of 5%. We have 
also briefly discussed the impact of our calculation on the 
experimental determination of the gaugino masses, but 
leave a detailed experimental simulation to future work. 
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